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We develop a search method for gravitational ringing of black holes. The gravitational ringing is due to complex 
frequency modes called the quasi-normal modes that are excited when a black hole geometry is perturbed. The 
detection of it will be a direct confirmation of the existence of a black hole. Assuming that the ringdown waves are 
dominated by the least-damped (fundamental) mode with least imaginary part, we consider matched filtering and 
develop an optimal method to search for the ringdown waves that have damped sinusoidal wave forms. 

When we use the matched filtering method, the data analysis with a lot of templates required. Here we have to 
ensure a proper match between the filter as a template and the real wave. It is necessary to keep the detection 
efficiency as high as possible under limited computational costs. 

First, we consider the white noise case for which the matched filtering can be studied analytically. We construct an 
efficient method for tiling the template space. Then, using a fitting curve of the TAMA300 DT7 noise spectrum, we 
numerically consider the case of colored noise. We find our tiling method developed for the white noise case is still 
valid even if the noise is colored. 



I. INTRODUCTION 



There are many on-going projects of gravitational wave detection in the world; LIGO[l], VIRGO [2], GEO- 
600J31, ACIGAM and TAMA300[5] which are ground-based laser interferometers, and EXPLORER[6l, ALUEGROEI, 
NIGBE[8], NAUTILUS [9] and AURIGAHO] which are bar detectors. Furthermore, there are some fixture space intfer- 
feromete? projects suchTIs LISAflll. 

Detection of gravitational wa\%5"provides us with not only a direct experimental test of general relativity but also 
a new window to observe our universe. To use them as a new tool of observation, it is necessary to know theoretical 
waveforms. Once we know them, we may appeal to the matched filtering technique to extract source's information 
from gravitational wave signals. However, because the signals are expected to be very weak and the amount of data 
will be enormous for long-term continuous observations, it is essentially important to develop efficient data analysis 
methods. 

For the ground-based and future space-based interferometers, the coalescences of compact object binaries are the 
most important sources of gravitational waves. The process of coalescence can be divided into three distinct phases. 
During an inspiral phase, the gravitational radiation reaction timescale is much longer than the orbital period. The 
gravitational waves from the inspiral carry the information of the masses and spins of the systems and so on. After the 
inspiral, compact object binaries encounter the dynamical instability and then would merge. This phase is called as a 
merger phase. The gravitational waves from the merger give us the information about the highly nonlinear dynamics 
of relativistic gravity. Finally, if a black hole is formed by merger, this system can be describe as oscillations of this 
final black hole's quasinormal modes and then settles down to a stationary Kerr state. The emitted gravitational 
waves in this ringdown phase carry the information about the mass and spin of the final black hole. 

In this paper, we consider gravitational ringing of distorted spinning (Kerr) black holes. The ringdown waves are due 
to quasi-normal modes of black holes that are complex frequency wave solutions of the perturbed Einstein equations 
with purely outgoing-wave boundary condition at infinity and ingoing-wave at horizon, with vanishing incoming-wave 
amplitude. A quasi-normal mode is characterized by the central frequency f c , usually called the (quasi-)normal-modc 
frequency, and the quality factor Q which is inversely proportional to the imaginary part of the complex frequency. 
The gravitational waves emitted at the last stage of the formation of a black hole are also expected to be dominated 
by the quasi-normal modes. 
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The quasi-normal modes can be obtained by solving the Teukolsky equation that governs the perturbation of a 
Kerr black hole. Their properties were analyzed extensively by Leaver [121, and it is known that the least-damped 
mode belongs to the I — m = 2 spin-2 spheroidal harmonic modes. The oependence of the parameters {/ c , Q} of the 
least-damped mode on the black hole parameters {M, J}, where M is the black hole mass and J is the spin angular 
momentum, is briefly reviewed in Sec. [H] Assuming that the ringdown waves are dominated by the least-damped 
mode, the black hole parameters {M, J} can be uniquely determined by measuring the parameters {/ c , Q} of the 
gravitational wave signal. 

It is noted that the ringdown signal decays exponentially. If the signal to noise ratio (SNR) is large, not much effort 
is needed to detect the signal, but if the SNR is small, the signal is going to be buried inside the noise even more 
after just one wave length, and there is absolutely no hope of seeing many more cycles, specifically because the signal 
drops exponentially with time. So, we have to treat the loss of the SNR not to be large. 



Methods to search for the ringdown waves were discussed previously by several authors. Echeverria [13] investigated 

the problem of extracting the black hole parameters from gravitational wave data in the case when the~5lgnal-to-noise 
ratio (SNR) is large. Finn [14] improved this situation by developing a maximum likelihood analysis method that 
can deal with any SNR. Flarragan and Hughes then considered the parameter extraction from the three stages of 
a binary coalescence, i.e., from inspiral, merger and ringdown phases, in their series papers [15] . For the ringdown 
phase, they discussed the relation between the energy spectrum of the radiation and the SNR. CTeighton [16] reported 
the results of analyzing data of the Caltech 40m by matched filtering, and emphasized the importance orroincidence 
event searches to discriminate spurious events from real events. But the search was limited to a single ringdown 
wave template. In order to treat ringdown waves with unknown parameters, we need to prepare a lot of theoretical 
templates. It is necessary to keep the loss of the SNR as small as possible under limited computational costs. So, we 
should consider an effective template spacing. Recently, Arnaud et al. [17] have discussed a tiling method to cover 
the 2-dimensional template space {/ c , Q}. In this paper, we develop a Tttfterent tiling method which is much more 
efficient than that of [171 and examine the efficiency with TAMA300 DT7 noise spectrum. 

Here we make a comment on combining this search technique of ringdown waves with the current search techniques 
for the earlier two stages of the process in the coalescence of compact object binaries, i.e., inspiral and merger phases. 
Ultimately, the search method and the data analysis approach should be able to handle all the three phases in a 
unified manner, going smoothly from one phase to the next. In the present state, it may be difficult to obtain the 
parameters of the final black hole from the gravitational waves in the inspiral and merger phases. However, in the 
case when a compact star is inspiralng into a super massive black hole, we will be able to obtain the information 
about the mass of a super massive black hole from the gravitational waves in the inspiral phase. This will improve 
the detection efficiency of ringdown waves significantly. 

The pape r is organized as follows. In Sec. [HI we briefly review the quasi-normal modes of spinning black holes. 
In Sec. IIIII we conside r th e template space {/ c , Q} for the white noise case analytically and develop an efficient 
tiling method. In Sec. II VI by using a fitting curve for the TAMA DT 7 noise spectrum , w e show that our tiling 
method developed in the case of white noise is valid even in the case of colored noise. Sec. is devoted to summary 
and discussion. In Appendix A, we discuss simpler cases when only cosine or sine part of the ringdown waves are 
considered. In Appendix B, we discuss the number of templates in the case when the mass of the black hole is known. 



A black hole is characterized by its mass M and spin angular momentum J. Here we use the dimensionless spin 
parameter a = J/Al 2 that takes a value in the range [0,1) with a = corresponding to a Schwarzschild black hole 
and a = 1 to extreme Kerr black hole. Quasi-normal modes of a black hole are complex frequency solutions of the 
Teukolsky wave equation that satisfy purely outgoing- wave boundary condition at infinity and ingoing- wave at horizon 
with vanishing incoming- wave amplitude. For fixed spheroidal harmonic indices (£,m), there are infinite number of 
quasi-normal modes. They are assigned with an index n with the order of the magnitude of the imaginary part, 
i.e., the n = 1 mode has the smallest imaginary part (or the largest quality factor). We call it the least-damped 
(fundamental) mode. It is known that the imaginary part of the I = m = 2 least-damped mode is the smallest of all 
the quasi-normal modes, and results of black hole perturbation calculations as well as numerical relativity simulations 
strongly suggest that the ringdown waves are dominated by this t = m = 2 least-damped mode unless the spin 
parameter a is extremely close to unity. Hence, we focus on this mode. 

For the I = m = 2 least-damped mode, analytical fitting formulas for the central frequency and quality factor for 
a black hole of mass M and dimensionless spin a were found by Echeverria as 




II. QUASI-NORMAL RINGING MODES 




(2.1) 
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2.0(1 - a) 
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(2.2) 



The ringdown waveform is expressed as 



Hfc Q, to, 4>o; t) = 



e -*fc(t-t )/Q cos (2 7r/ c (t 




*o) - <Ao) for t > t 
for t < to 



(2.3) 
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where to and fa are the initial time and phase of the ringdown wave, respectively. 



III. TEMPLATE SPACE 



In this section, we develop an efficient technique of template spacing which can be used for matched filtering of the 
quasi-normal ringing waveforms. Here, the detector noise is assumed to be white noise to make it possible to deal 
with the problem analytically. The case of the colored noise is discussed in the next section. 



A. Distance function 

We have temporarily set the amplitude to unity for simplicity in Eq. (|2.3|) . Note that the knowledge of the amplitude 
is not necessary for the template spacing in matched filtering. For t>to, the ringdown wave (|2.3|) is divided into two 
parts. 

Kfc, Q, h, fa; t) = h c (f c , Q, t ; t)cosfa + h s (f c , Q, t ; t) sin fa , (3.1) 

where 

M/« Q, to; t) = e -" /c v~ t0) cos(27r/ c (t-t )), (3.2) 
h s (f c ,Q,t ;t) = e -" fc Q~ tn) sin(27r/ c (t-t )). (3.3) 

Performing the Fourier transformation h c / s (f) = j a ° oo dte 2 ' K ' t * t h c / s {t), we obtain the waveform in the frequency 
domain as 

~ h(f Off)- (fc-2ifQ)Qe* i *f t ° 

cU ° y ' " n(2f c Q-if c -2fQ)(2f c Q + if c + 2fQ)> [6A) 

2 f c q 2 e 2 i f t 

hs{f c , Q, to; f) = w{2fcQ _ ifc _2 fQ){2fcQ + ifc + 2fQ) ■ ( 3 - 5 ) 
The waveform in the time domain is real, so the following relation is satisfied. 

h*(f) = h(-f), (3.6) 

where the star ( *) denotes the complex conjugation. 

Here, we introduce the inner product between two functions as 

(a, 6)= J ™" dfa(f)b*(f), (3.7) 

where / ma x is the maximum frequency we take into account in the analysis. In the actual data analysis, it is equal to 
or less than the half of the sampling frequency of data. 

In the matched filtering, we calculate the inner product between the template h and the signal x defined by (x, h). 
We first introduce the normalized template. We define the normalization constants as 



N c {fc, Q, t ) = (h c (fc, Q, t ),h c (fc, Q, t a )) 

1 (2Q 2 + 1)Q 

2 tt(4Q 2 + 1)/ c ' 

N s (f c , Q, t ) = (h s (f c ,Q,t )Mfc,Q,to)) 
tt(4Q 2 + 1)/ c ' 

when we set / max = oo. The normalized templates h c / s (f c , Q, to; f) are given by 



(3.8) 



(3.9) 



hc/ s {fc, Q, h; f) = ==h c/s (f c , Q, to; f) . (3.10) 

v iV c/sUc y, to) 



4 



We then consider h = h c cos <po + h s sin^o as a template. We note that the h c and h s are not orthogonal. Their inner 
product is given by 

(h c (fc, Q, *o)A(/c> Q, t )) = 



V 2 (2Q 2 + 1) 
=: c(/c, Q, to), (3.11) 

when / max = 00. In this case, the maximization of (x,h) over the phase 4>o can De carried out analytically to yield 

m 

A(ir) = max(:c, h) 

(x, h c (f c , Q, t )) 2 + (x, hjjc, 0, to)) 2 - 2c(/ c , Q, t Q )(x, h c (f c , Q, tp))(x, h s {f c , Q, t )) 

l-c(f c ,Q,t ) 2 ' 

In the following, we consider the 3-dimcnsional template space {/ c , Q, to}. 

Here we consider the match C(df c , dQ, dto) between the template with (f c , Q, to) and the normalized signal having 
slightly different sets of the parameters (f c + df c , Q + dQ, to + dto). Since we have already maximized over the phase 
<f>o in A, without loosing generality, we only need to consider h c for the signal. Then, the match is defined by 

C(df c , dQ, dt ) = A(h c (f c + df c , Q + dQ, t + dt )) 

= 1 2 ^ 4 2Q 2 (4Q 2 + 5) 2 

(2Q 2 + l)/2 Jc (4Q 2 + 1) 2 (2 Q2 + 1) H 

2Q 2 , n 

/ c (4Q2 + i)( 2 Q2 + 1) ^^ 

, 4tt/ c (1 + 4Q 2 )/ max 2 tt/ c 

+ Q(2Q2 + 1) *o + 2g2~j7Y d<0 ^ ' ( 3 - 13 ) 

where we have used the approximation / max 3> / C) and only the leading term in / max //c are shown. We will find that 
the dependence of / max will be cancelled out and does not appear in the final result. The inequality C(df c , dQ, dt ) < 1 
means that there will be a loss of signal-to-noise ratio unless the actual parameters of a gravitational wave signal fall 
exactly onto one of the templates. 

The smaller the match C is, the larger the distance is between the two signals in the template space. Therefore, 
we define the metric in the template space by ds 2 ^ = 1 — C [19], that is, 



dsf 3 \ = g^dx i dx :> 



'(3) - y« 



2Q4 



(2Q2 + l)/2 (4Q2 + 1)2 (2 Q2 + 1) 

2Q ' df c dQ 



/ c (4Q2 + i)( 2 Q2 + i) 
, 47r/ c (l + 4Q 2 )/ max 2 irf c 

+ Q(2Q2 + 1) dt ° + dt ° dQ ■ (3 ' 14) 

B. Projection to 2-dimensional template space 

Now to maximize the match with respect to dto, we consider a projection of the distance function into two dimensions 

(3) 

spanned by f c and Q. Namely, we project g\-' to a two dimensional subspace orthogonal to the to-axis ^| as 

(3) (3) 

(2) (3) 9n 9jt ,„..-,. 

9lJ — 9l,J (3) ) (O.IOJ 

9t t 



where the indices /, J are {f c , Q}. So, we find 

5 (2) 



ds 2 2 \ = gj]dx I dx J 
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2 Q 4 „ a , 2Q 2 (4Q 2 + 5) 2 



(2Q 2 + 1)/| Jc (4Q 2 + 1) 2 (2 g 2 + l) 
2Q 3 



/ c (4Q2 + i)(2Q2 + i; 



d/cdQ, (3.16) 



where we have taken the limit / max — ► oo. 

It is noted that the cross term <i/ c ciQ arises here. Later, we discuss a coordinate transformation in the template 
space that removes the cross term, in order to make our analysis of the template spacing and the error estimation 
easier. Wc also note that the dependence of the metric on f c can be eliminated by the simple coordinate transformation 
fc -> ln/ c . 

It is also noted that the smaller the volume element of the metric is the fewer the number of required filters is to 
cover the template space. 

The metric in the case of only h c or h s as templates is discussed in Appendix 1X1 



C. Diagonalization of the metric 

Let us perform the coordinate transformation by which the two dimensional metric (|3.16|) is transformed to a 
diagonal, conformally flat metric. 

We start from the coordinate transformation that removes the frequency dependence in the metric. Wc set 



dF=^, (3.17) 

Jc 



which gives 



ds z = g FF dF z + g QQ dQ z + 2g FQ dFdQ 



(2Q 2 + 1) (4Q2 + 1)2( 2 Q2 + 1 ) 

- W+ %Q* + l) dFdQ - (3 ' 18) 

The transformation that removes the off-diagonal element is found by setting F = X — u(Q) and requiring 

9ffu'(Q)- g FQ = 0. (3.19) 

We find 

HK 4+ ^)' (3 ' 20) 

which gives 

ds 2 = g FF dX 2 + (g FF u' 2 - 2g FQ v! + g QQ ) dQ 2 

= g FF ( dX 2 + 9FF9QQ 2 - g ^ dQ 2 \ . (3.21) 
\ 9ff J 

Then we can perform a further coordinate transformation to make the metric conformally flat. Namely, by the 
transformation Q — > Y defined by 



Jq V 9ff{Q') Jq2 x(4x+1) 

we obtain 

ds 2 = Q(Y) (dX 2 + dY 2 ) , (3.23) 

where the conformal factor is given by 

Q(Y)=g FF (Q(Y)). (3.24) 
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Here Q is now a function of Y determined by inverting Eq. (|3.22|) . 

Although the above coordinate transformation involves complicated functions that may not be expressed in terms of 
elementary functions, we find it is sufficient to use their large Q expansion forms for Q > 2. Up to 0(1/Q a ) inclusive, 
we have 



11 1111 11 
8Q 2 ~ 64 Q4 + qg ~ 2048 Q« 
I I 11 3 1 11 17 1 



X = F + u = F + ln2 + -^--- i + — (3.25) 



Y = I I (3 26) 

2Q 24 Q 3 160 Q 5 128 Q 7 4608 Q 9 ' y ' ' 



To the same accuracy, the inverse transformation becomes 



Q= II + i F _5! F 3 + 166 yt _59I7 F r. (328) 

^ 2 Y 6 90 135 1350 v ; 

The conformal factor fl(Y) is given by 

il{Y} - 4i^~3 + 60 y "54 y + Yl00 Y ■ (3 - 29) 

When Q = 2, the errors induced by the above expansion are found to be ~ 0.1 %. This is accurate enough for our 
purpose as long as we allow the SNR loss, ds^ ax , of a few percent. 



D. Tiling method: basis 

In the previous subsection, we have derived the simple, conformally flat metric (|3.23J) for the template space. Here, 
using this metric, we formulate a tiling algorithm which is not only efficient but also quite simple. 

To develop such a method, we note the following. Because of the conformal flatness, the contour of the fixed 
maximum distance ds 2 = ds^ax centered at a point on the (X, Y")-plane is a circle for sufficiently small ds^ ax . 
Furthermore, along a line of Y =constant, Q(Y) is constant. Thus, choosing first an appropriate Y =constant line, 
say Y = q±, we may place circles of the same radius with their centers located along the line Y = q\ to cover a region 
surrounding that line. Then, if we find an algorithm to place circles along the Y = q± line and another algorithm 
to choose the next Y =constant line, say Y = q2, to be covered in an appropriate way, we can repeat this tiling 
procedure to cover the whole template space. 

Let us assume that the template space to be tiled is a rectangle given by F m { n < F < F max and Q m ; n < Q < <5 max . 
In the (X, Y) coordinates, this rectangle is mapped to the region bounded by the two Y ^constant lines corresponding 
to Q = Qmin and Q = Q max , which we denote by Y = Yq and Ym, respectively, and the two lines X = v m - m (Y) and 
X = u max (y) corresponding to F = F min and F = F max , respectively. Note that Y > Ym since large Y corresponds 
to small Q. 

First, we construct a method to determine the spacing of the circles along each Y =constant line. Let us consider 
the line Y = q and place two circles with radius r centered at (p, q) and (p + Ap, q), 

(X -pf + (Y -qf =r 2 , (X -p- Ap) 2 + {Y -q) 2 =r 2 . (3.30) 

We assume Ap < 2r so that the two circles intersect at the two points, (p + Ap/2 , q ± d(r; p)), where d(r; p) is the 
distance to each intersecting point from the line Y = q (see Fig^l, given by 



d(r; Ap) = ^r 2 -^-. (3.31) 

Our purpose is to tile the template space by the smallest possible number of filters. In order to do so, we choose the 
parameter p in such a way that the area defined by S — Apd(r; Ap) is maximized, i.e., 

Ap = V2r, (3.32) 

d = . (3.33) 

The radius r is determined by the value of ds^ M and q as 

r2 = d -m- (3 - 34) 
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In this way, we tile the region that covers the line Y = q. 

To choose the first line to be covered, we start from the point on the (X, Y)-plane corresponding to (F, Q) = 
(F max , Qmin), that is, (X, Y) = (Xq, Yq) where X = v max (Y ). Then we choose the first line Y = q\ and the radius 
7*1 so that the point {Xq,Yq) is just on the edge of the first circle and an intersecting point of the first and second 
circles lies on the line Y = Yq as Fig. [21 This is achieved if the center of the first circle is at 

(3.35) 



(3.36) 



(pi, Qi) = [Xq - n/V2, Y - n 



with the radius determined by solving Eq. 



dsl 



, which reads in the present case, 
= r\ O(Y - n/V2) . 



(h, YJ 



0.25 



0.15 




0.05 



FIG. 2: Choosing the first line on the (X , Y)-plane. 



Then the center of the n-th circle is at 

(p n , qi) - (x Q - (2n - 1) n/V2, Y - n/V2j 

and the number of circles needed to cover the line Y — q± is given by 

Li 



V2 



7*1 



(3.37) 



(3.38) 
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where [a;] denotes the maximum integer smaller than x + 1 and L\ is the coordinate length of X to be covered, i.e., 



Li = v mttX (Y ) - v min (Y ) 



(3.39) 



Once the covering of the first line is done, the second Y =constant line is chosen as follows. Let Y\ = Yq — V2r\ 
and let (Xi,Yi) be the intersecting point of the lines Y — Y\ and X — v max (Y), i.e., (Xi,Yi) = (v max (Yi),Yi). We 
choose this point as the starting point for the covering of the second line as Fig. [21 That is, the center of the first 



0.25 



0.15 




0.05- 



FIG. 3: Covering of the second line. 

circle {p 2 ,<l2) on the second line Y = q 2 is 

(pa, ©) = (Xi - r 2 /V2, Yi - r 2 /V2) , (3.40) 

with the radius r 2 determined again by Eq. (|3.34|) . Then the second line is covered by the same procedure we took 
for the first line. We repeat this procedure until we tile the whole template space we need to cover^ 

With this tiling procedure, the total number of templates is given as follows. Generalizing Eq. (|3.38|l . the number 
of templates for the i-th Y =constant line (Y — g,) is 



Ni = 



L 



y/2l 



where is determined by 



dst aax = rfn(Y i ^-r i /V2 



(3.41) 



(3.42) 



The number of Y =constant lines necessary to cover the template space is determined by the minimum integer v that 
satisfies the inequality, 



Y,V2n>Y -Y M . 



(3.43) 



And the total number of templates is 



(3.44) 



9 



E. Tiling method: application 

Let us apply the method developed in the previous subsection to the case of the parameter space (f c , Q) which has 
the range, 

10 2 Hz < f c < 10 4 Hz , 

2 < Q < 20. (3.45) 

We set cLs 2 lax = 0.02. This choice is made in order to make the SNR loss to be ~ 3 % in the presence of colored noise 
as discussed in the next section. Using Ea. H3.17|) . we set F = ln(/ c /100Hz). The above range corresponds to 

6.93 x 10" 1 < X < 5.33, 

2.50xl0~ 2 < Y <2.55xl0 _1 . (3.46) 

in the (X, Y~)-space. 



(fc Q) 


(X, Y) 


10 2 Hz,2.0 


7.234594008 10" 1 , 2.546762255 10" 1 


10 2 Hz, 20.0 


6.934595830 10 _1 , 2.500520248 10~ 2 


10 3 Hz,2.0 


3.026044494, 2.546762255 10 _1 


10 3 Hz,20.0 


2.996044676, 2.500520248 10~ 2 


10 4 Hz,2.0 


5.328629588, 2.546762255 10" 1 


10 4 Hz, 20.0 


5.298629769, 2.500520248 10~ 2 



TABLE I: The relation between the two coordinates 

Following the procedure described in the previous subsection, we choose the starting point (Xq, Yq) which corre- 
sponds to (f c , Q) = (10 4 Hz,2.0). The radius r\ is determined by Eq. (|3.36l) . or 

0.02 = r 2 fi(y -r 1 /v / 2). (3.47) 
Then all the parameters for our tiling method are determined. In Table IITI we summarize the tiling parameters. 



i 


(Pi, ?i) 


(Xi, Yi) 


Ti 


Ni 







5.328629588, 2.546762255 10" 1 






1 


5.285172504, 2.112191414 10" 1 


5.311957476, 1.677620573 10" 1 


6.145759769 10" 2 


53 


2 


5.283700230, 1.395048115 10" 1 


5.304418737, 1.112475658 10" 1 


3.996178019 10 -2 


82 


3 


5.285789828, 9.261865694 10" 2 


5.301037366, 7.398974808 10" 2 


2.634525556 10" 2 


124 


4 


5.288679812, 6.163219420 10" 2 


5.299527944, 4.927464032 10~ 2 


1.747622029 10" 2 


187 


5 


5.291307826, 4.105452276 10" 2 


5.298855740, 3.283440519 10~ 2 


1.162500174 10" 2 


281 


6 


5.293381065, 2.735973019 10" 2 




7.742359638 10" 3 


421 



TABLE II: The parameters for the tiling of the template space are summarized. {Xi, Yi) denotes the starting point of the i-th 
template spacing along the line y = qi. (pi, qi) denotes the center of the first circle along the i-th line y — qi, rt is the radius 
of the circle, and Ni is the number of circles necessary to cover the i-th line. 

It is noted that the number v of the Y ^constant lines is very small, 



v = 6 



(3.48) 



10 




' 5.15 52 ' ' 5.25 ' ' ' 5.35 

X 



FIG. 4: A part of the tiling of the template space in the (X, Y) coordinates. Templates are taken at the centers of the circles. 

0.3" 
0.25" 

0:.2- 

Y 0.15- 

0.1 - 



0.05 



— , , , , , , 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ 

J 1 2 3 4 5 

X 

FIG. 5: The template points on the (X, F)-plane. It is noted that we may put the templates with the same interval along the 
X direction. 



The total number of templates is calculated to be Af = 1148 (J\f — 1780 for Q niax — 34.3). In Fig. 21 we show the 
tiling of the template space in the (X, Y) coordinates. The tiling of the template space in the original coordinates 
(/cj Q) is shown in Fig.EJ 

In order to see how effective this tiling method is, we calculate the ratio between the sum of the areas of all the 
circles <S c ir and the area to be covered Spar- We find 

*^cir 

11 = s~ 

'-'par 

= 1.57. (3.49) 

Here we have adopted Q m ax = 22.9 which is slightly larger than the pre-assigned value of Q max = 20. The former 
value corresponds to the value Y — qe — tq/\/2, i.e., the value of the Y =constant line marginally covered by the 
circles centered along the line Y = qg. This ratio corresponds to Tftot m the paper by Arnaud et al. H7], in which they 
obtained Tjtot = 2.12 for their tiling method. Thus our tiling method is far more efficient (and simpler) than that of 
Arnaud et al.. 



IV. VALIDITY OF TILING METHOD IN THE CASE OF COLORED NOISE SPECTRUM 

The tiling method discussed in the previous section is based on an assumption that the noise is white. We expect 
this assumption is very good because the ringing wave is rather narrow banded except the case Q ~ 2. In order to 
confirm this, we examine the effectiveness of the tiling method in the case of colored noise. 
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FIG. 6: A part of the tiling of the template space in the original coordinates (/ c , Q) with f c measured in units of 100Hz. Note 
that the contour of ds^ ax = 0.02 for each template is warped. 
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FIG. 7: The template points on the original coordinates (f c , Q). The smaller the value of / c is, the larger the number of 
templates is needed. 



As a model of detector's noise, we use a fitting curve of the one sided noise power spectrum of TAMA300, which is 
given by 

* n** f &5 \ 63 1 /220\ 10 1 /710\ 3 3 If f \ 2 1 / / \ 6 

This formula of the noise spectrum is obtained by fitting Fig. 3.4.1 in TAMA REPORT 2002 [20] which is based on 
the spectrum during Data Taking 7 in 2002, and is valid between 60 Hz and 40000 Hz. Here~We have ignored the 
overall amplitude of S n , because it does not affect the results. 

We prepare the template bank using the analytical method in the white noise case. The minimal match is assumed 
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to be 0.98. We also generate signals which are normalized to unity. We then calculate the maximum of the match 
between the signal and templates. When the signal is completely the same, the match becomes unity. However, we 
have the match less than unity due to mismatch between signal and templates. If the match is always greater than 
0.98, we have a justification to use the tiling method even in the case of the colored noise spectrum. 

We used 2500 signals. The range of the central frequency f c and the Q value of the signals are 100Hz < f c < 10 4 Hz 
and 2 < Q < 20 respectively. The value of f c and Q are randomly given in this region. 

In Fig. |8l we show the number of signals in terms of the value of the match. 

3501 I 1 



300- 



250: 



200- 



150: 



100- 



50: 



0.98 0.985 0.99 0.995 1 

FIG. 8: The number of signals in terms of the value of the match between the signal and templates. The bank of templates 
are determined assuming the minimum match 0.98 using the method in Section III. The mean value of the match is 0.993. 

From Fig. [SJ we find that we can detect the most of the signals without losing 2% of the signal-to-noise ratio by 
the tiling method assuming white noise. Thus, the tiling method, constructed analytically assuming white noise, is 
valid even in the case of the TAMA noise spectrum. 



V. DISCUSSION 

In this paper, we proposed a tiling method of the template space in the matched filtering search for the gravitational 
waves of black hole ringing. 

First, we discussed a tiling method assuming the detector noise is white. We found that the l/Q-expansion was 
useful in performing the coordinate transformation of the metric of the template space because of Q > 2. We also 
find an efficient tiling method of the template space which can been formulated analytically. When the range of f c 
and Q is 10 2 < f c < 10 4 Hz and 2 < Q < 20, and the minimal match is 0.97, the ratio rj, between the sum of the 
areas of the equal match circles around the templates, 5 c i r , and the area to be covered, <S par , is 1.57 which is much 
smaller than the value obtained before. 

Next, we have discuss the validity of this tiling method in the case of colored noise spectrum. As a model of realistic 
noise power spectrum, we used a fitting curve of the noise power spectrum of TAMA300 during 2002. We prepare a 
template bank assuming the minimum match 98%. We then examined the loss of signal-to-noise ratio of the signals 
detected in the template bank. We found that, in most case, we do not lose the signal-to-noise ratio no more than 
2% which is expected from the pre-assigned minimum match. This shows that the tiling method can be used even in 
the case of colored noise spectrum. 

The tiling method should be tested using the real data of laser interferometers. A study in this direction using the 
data of TAMA300 is now progressing [211 . 

In the analysis using real interferometers' data, we have to treat the non-stationary, non-Gaussian noise. It is found 
[221 that we observe even more fake events than in the case of inspiraling wave, since the ringing wave is usually much 
shorter than the inspiral waves and can easily be affected by short bursts. In this situation, we would need to introduce 
some methods to remove such fake events without losing real ringing gravitational wave signals. Coincidence analysis 
between several detectors would also be needed to reduce the fake event rate. Further, since the ringing waves may 
be excited by the detector itself, we may need to perform coincidence analysis to reject such spurious events anyway. 
We will also work on this problem in the future. 
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APPENDIX A: THE METRIC IN SOME SPECIAL CASES 



In this appendix, we summarize the distance function in the case when we consider only the cosine or sine part of 
the wave, ignoring the contribution of the phase and the initial time. We also consider the case when the detector 
noise is white. Although the applicability of such special cases will be limited, we describe them here because there 
are still some cases when these formulas become useful. 



as 



1. Cosine wave case 

We consider the cosine waveform ignoring the initial time and phase of a ringdown wave; 

h(/c) g. t) = (e-*A*/« cob(2 */„*) fort>0 , (A1) 
( lor t < . 

Performing the Fourier transformation h(f) = dte^^h^), we obtain the waveform in the frequency domain 

~ hif Q.f)- (/c-2»/Q)Q (A2) 

In the following, we consider normalized waveforms. The normalization constant is derived as 

N(f c , Q) = ^ \h(f c , Q; f)\ 2 df 

J — oo 



_ 1 (2Q 2 + 1)Q 
2 ^(4Q 2 + 1)/ C • 

With the above normalization constant, the normalized wave h(f c , Q; /) is given by 



(A3) 



h(f c , Q; /) = —=L=h(f c , Q; /) . (A4) 

Now, we consider the correlation between the two normalized waves having slightly different sets of the parameters 
(f c , Q) and (/c + rf/c, Q + dQ), 

C(df c , dQ) = / df h(f c , Q; f)h*(f c + df C7 Q + dQ; f) 



_ 1 16Q 4 + 6Q 2 + 1 2 _ 1 64Q 8 + 128 Q 6 + 28Q 4 + 1 2 
8 (2Q 2 + 1)/ C 2 /c 8 Q 2 (4Q 2 + 1) 2 (2Q 2 + 1) 2 Q 
1 (8Q 4 + 2Q 2 + 1) 
+ 4/ c Q(4Q 2 + l)(2Q 2 + l) dfcdQ - (A5) 

The inequality C(df c , dQ) < 1 means that there will be a loss of the signal to noise ratio unless the actual parameters 
of a gravitational wave signal fall exactly onto one of the templates. 

The smaller the correlation C is, the larger the distance is between the two signals in the template space. Therefore, 
we define the metric in the template space by ds 2 = 1 — C, that is, 

2 1 16Q 4 + 6Q 2 + 1 2 , 1 64Q 8 + 128 Q 6 + 28Q 4 + 1 2 
S ~ 8 (2Q 2 + 1)/ 2 /c + 8 Q 2 (4Q 2 + 1) 2 (2Q 2 + 1) 2 Q 

1 8 « 4 + 2 « 2 + 1 dfcdQ , ( A6) 



4 / c Q(4Q 2 + l)(2g 2 + l) 
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Note that the dependence of the metric on f c can be eliminated by the simple coordinat e tr ansformation f c — > ln/ c . 
The metric in the case of the sine-wave templates instead of the cosine-wave templates (|A1J| is calculate in the next 
subsection. For large Q, the two metrics coincide with each other, implying that the phase effect is weak in the large 
Q limit. 

The number of filters, Nf, which are required to achieve the detection efficiency determined by the maximum 
allowable distance ds max (i.e., the maximum allowable loss of SNR) in the region / m in < / < /max and Q m i n < Q < 
Qmax is estimated by integrating the volume element y/det g over that region of the template space. We find 



N f = T^max) 1 ln(/max//min) 



dQ 



Q(16Q 6 + 32 Q 4 + 10Q 2 + 1) 



1/2 



(2Q 2 + l) 3 / 2 (4 Q 2 + l)!/2 

The Q m ax-dependence of Nf is shown in Fig. For large Qmax, it can be derived analytically as ^ 

V2, 



N 



f I large- Q 



: ( ds max) 1 ln(/max//min) 1 



(A8) 



(A9) 



6 8 10 12 14 16 18 20 



FIG. 9: The Qmax-dependence of the number of required filters [2ds2iax-^f/log(/max//mm)]. Here we chose Q m i n = 2. 

From Fig. |51 we see that the large Q nl „ y approximation is good even for smaller Q max as far as the number of 
required filters is concerned. In Table ITTTT we show the number of filters in the case Q m i n = 2, / m ; n = 100Hz, and 
^ s max = 0-03 for several values of / max and Q ma x- We also note that the number of required filters depends on the 
choice of a tiling method in practice. 



/max (Hz) 


Qmax — 10 


Qmax — 10 


Qmax — 10 


10 3 


218 


2,661 


27, 084 


10 4 


436 


5,321 


54, 167 



TABLE III: The estimation of the number of filters with the maximum SNR loss of 3 % 



Next, let us perform the coordinate transformation by which the cosine-wave template metric (|A6|) is transformed 
to a diagonal, conformally flat metric. 

We start from the coordinate transformation that removes the frequency dependence in the metric, 



which gives 



dF = d\nf c , 
ds 2 = g FF dF 2 +g QQ dQ 2 + 2g FQ dFdQ 



(A10) 
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1 16Q 4 + 6Q 2 + 1 a 1 64 Q 8 + 128 Q 6 + 28 Q 4 + 1 a 
8 (2Q 2 + 1) + 8 Q 2 (4Q 2 + 1) 2 (2Q 2 + 1) 2 g 

1 8Q 4 + 2Q 2 + 1 

"4 Q(4Q 2 + l)(2Q 2 + l) ^ dQ - (AU) 
The transformation that removes the off-diagonal element is found by setting f = X — u(Q) and requiring 

g FF u'{Q)-g FQ = 0. (A12) 

We find 



u 



3y/7 
14 

which gives 



arctan 



3 + 16Q 2 \ 7T 



\/7 / 2 



1 16Q 4 (16Q 4 + 6Q 2 + 1) 
~4 ln ' (A13) 



ds 2 = «? FF dX 2 + (,g FF u' 2 - 2.g FQ u' + g QQ ) dQ 2 = g FF I dX 2 + 9ff 9Q ® 9f q \ d Q 2 . ( A 14) 

\ 9 FF J 

Then we can perform a further coordinate transformation to make the metric conformally flat. Namely, by the 
transformation Q — > Y defined by 



v r ^ det 9(Q') r a 2V64a; 4 + 144x3 + 72x 2 + 14^+1 

we obtain 

ds 2 = 0(7") (d^ 2 + dF 2 ) , (A16) 

where the conformal factor is given by 

Sl(Y)=g FF (Q(Y)). (A17) 

Here Q is now a function of Y determined by inverting Eq. (|A15|I . 

Although the above coordinate transformation involves complicated functions that may not be expressed in terms of 
elementary functions, we find it is sufficient to use their large Q expansion forms for Q > 2. Up to 0(1/Q S ) inclusive, 
we have 

_L__L_ _3 1_ _L3__|_ 43 1 109 1 

+ U ~ + 16 Q 2 ~ 256 Q 4 + 3072 Q6 ~ 32768 Q 8 + 327680 Q 10 ' ( ' 

_ y/2 / 1 1 1 73 1 21 1 6127 1 \ 

~4~ \Q + 12Q3 ~ 640Q5 + 256 Qi ~ 98304 QV ' ^ ' 

To the same accuracy, the inverse transformation becomes 

f = x — -Y 2 + —Y i - ^ + ™r 8 - H^iVo (A20) 

2 12 45 360 2025 

V2 / 4401959 , 9892 . 737 , 



Q = — Y' + Y° Y A + -Y H . (A21) 

^ 4 V 5400 135 90 3 ry v ; 

The conformal factor O(F) is given by 

= _749639 y6 1255 y4 _ 119 J_ 1 1 
v 7 5400 108 120 24 8 Y 2 y ' 

The structure of this metric is the same as the one in the general case discussed in Section III. The effective tiling 
method describe in Section III can thus be used to set the template space using this metric. 
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2. Sine wave 

Next, we consider the damped sine wave e " sin(2 7r / c i) . The calculation is done as same as the cosine wave 
case. First, we obtain 

*" - i (8 (f + 3) + i g, 1 ^), «■ - J q^TJ « ■ (A23) 



For the above metric, the l/Q-expansion is performed and we obtain 

„2 _ „ f _ jt?2 i jr>2 

l2 



ds 2 = g FF dF 2 +g QQ dQ 2 + 2g FQ dFdQ : (A24) 



11 11 3 1 5 1 7 1 
8Q 2 " ~ 16 Q 4 + ge - 256 Q 8 + 1024 ' 
1111 11 11 11 
9FQ ~ ~8Q ~ T6Q3 + 64Q5~ 256 + 1024Q9 ' 

The expansion is considered up to 0(1 /Q s ) beyond the leading term. 
Next, we consider the following coordinate transformation, 

X = F — — — — 11 1 49 1 _ 179 1 

+ 16 Q 2 + 256 Q 4 3072 Q 6 + 32768 Q 8 327680 Q 10 ' ^ ' 



_ y/2 ( 1 1 1 27 1 43 1 1067 1 

~ ~ \Q~ 6 Q3 + 640 ~ 3584Q 7 + 294912 ' ' ^ ^ 



The inverse transformation of this is given by 

181651 









15263 


X - 


V 19 Y 4 


136 y6 




2 12 


45 


2520 


V2 


( 71033 r7 


1073 y5 




4 


V 37800 


945 


90 



14175 » v ; 

+ (A28) 

Finally, we obtain the conformally flat metric as 

ds 2 = ft(Y) (dX 2 + dY 2 ) , (A29) 

"^-iao^+ss^ + Bo^ + H + i^- (A30 » 

The structure of this metric is the same as the one in the the general case of Section III and the case of cosine wave. 

Especially, the leading term with respect to 1/Q are the same as in the case of cosine wave. Therefore, the effective 
tiling method described in Section III is also applicable to this case. 



APPENDIX B: THE NUMBER OF TEMPLATES WHEN THE MASS IS KNOWN 

In this appendix, we consider the case when a compact star with mass 1 ~ 10 2 M Q is inspiraling into a super 
massive black hole with mass 10 6 ~ 1O 8 M , and after the final plunge, the ringing wave is excited. In such cases, 
it is expected that masses of the super massive black hole and the compact star is determined accurately during the 
inspiral phase. We do not know exactly the mass of the final black hole after the plunge of the compact star because 
we do not know how much energy are radiated as the gravitational waves and how much mass loss will be occurred 
in the case of neutron star. In this case, however, we can assume the mass of the final black hole within the accuracy 
10 -8 ~ 10~ 4 since the mass ratio is very large. In this case, we can perform an analysis to search for the ringdown 
wave using these mass parameters. Here we investigate how much templates we need in this situation. 

When we know the mass M of a black hole, we only need to investigate the spin of the black hole. From Eqs. I|2.1|) 
and l|2.2|) . the relation between f c and Q becomes 



320 



1 - 0.63 



-2/3 



M 

«0 



(Bl) 
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where we normalize the frequency by 100Hz. Now, we consider the white noise case. We only consider the cosine 
wave case. These assumptions are sufficient for the purpose here. From Eq. (|B1(I . 

B f° 210 f— VV 5/3 . (B2) 



dQ \M Q 
The metric (|A6(1 in the large Q limit is derived as 



ds 2 ~ (0.45 Q" 4/3 - 0.17 Q- 8/3 + 0.13 Q~ 2 ) dQ 2 

~ 0.45 Q- i/3 dQ 2 . (B3) 

When the maximum loss of the signal to noise ratio ds 2 lax is given, the number of templates needed to achieve this 

efficiency is estimated as N ~ (Qmax — Q^ n ) / ds max . If we consider Q m - m = 2 and Q max = 100, we only need 20 

templates. Thus, the prior knowledge of the black hole mass will make the matched filtering analysis substantially 
easy. 
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